Polymer Stretching by Turbulence 
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The stretching of a polymer chain by a large scale chaotic flow is considered. The steady state 
which emerges as a balance of the turbulent stretching and anharmonic resistance of the chain is 
quantitatively described, i.e. the dependency on the flow parameters (Lyapunov exponent statistics) 
and the chain characteristics (the number of beads and the inter-bead elastic potential) is made 
explicit. 
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A great task is to develop a theory of turbulence for 
a dilute polymer solution which explains, particularly, 
polymer drag reduction, the so-called Toms' effect [Q. 
Although a very plausible qualitative theory of polymer's 
drag reduction, describing the bulk viscosity and the size 
of the viscous boundary layer increase was given by Lum- 
ley no quantitative theory of the phenomenon is 

available. It is also widely accepted that the enormous 
increase of the polymer contribution into the stress ten- 
sor (by factor of ten thousand in a turbulent flow, e.g. 
pj, and up to 130 times in the recent experiment of 
Groisman and Steinberg on elastic turbulence at very 
low Reynolds number ||)) is associated with extreme 
polymer stretching J|,|| (there exists an alternative the- 
ory suggesting that the connection between the turbulent 
strain rates and the elongational viscosity is circumstan- 
, see also a comprehensive discussion of the subject 
|). Therefore, it seems logical to start a theoretical 
exploration of the whole problem of the polymer solu- 
tion hydrodynamics from description of a single polymer 
stretching. 

In the present letter I discuss a polymer chain placed in 
a statistically known random flow under the assumption 
that the back reaction of the polymer on the flow is neg- 
ligible, i.e. the polymer is passive. Polymer stretching 
in laminar flows was extensively studied in the past, e.g. 
P). A recent experimental advance has come with the 
state-of-the-art capability to visualize a polymer (DNA 
molecule) and therefore to study its stretching directly 
HI). A qualitative theory of polymer stretching by 
chaotic flow || has predicted the jump of the mean 
square polymer length at O(l) value of the stretching- 
to-linear elasticity ratio. Based on the recent progress 
in the theory of the passive scalar advection by a large- 
scale flow jlO] p"9| Balkovsky et al. [M have developed 
the approach of |3j] constructed a quantitative theory of 
the transition onset. In the present letter I generalize 
the approach of |^(|. Nonlinear, but still passive theory 
of polymer stretching above as well as below the afore- 
mentioned threshold is delivered . Of special note is the 



described dependence on the equilibrium polymer con- 
tour length. 

Equilibrium and even extended-by-flow lengths of a 
polymer chain are typically much smaller than the Kol- 
mogorov viscous scale. This allows to model a veloc- 
ity field (counted from a reference point, say one of 
the polymer's end points) advecting the polymer by 
u a (t;p) — a al3 (t)pP (i.e. the velocity at all scales un- 
der consideration changes linearly with the separation). 
Here the statistics of the traceless a (we consider incom- 
pressible flow) is assumed to be known. I will analyze a 
sequence of models moving steadily from ([!]), introduced 
in pCfl , through (Q) to (Jilf). Polymers will be modeled 
by a chain of beads (two beads in the cases of (|]) and 
(H) and N beads for ([l6])) connected by elastic (harmonic 
for (0) and anharmonic for @ and(|l6|)) springs. Forces 
of friction acting on a bead, which occurs as the bead is 
moving against the local flow, are equilibrated by elastic 
tension within the springs attached to the bead and by 
thermal forces. The goal is to describe statistics of the 
passive polymer(s) advected by u(t;r). 

I start by elaborating on the features and defects of the 
two-bead-one-linear-spring model. The basic equation of 
motion for the inter-bead separation p is 



dt 
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T a/3 



-kp a +£, a (t), 
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where the friction force on the rhs of (|l|) is equated by 
elastic and thermal forces on the lhs. In (|l|) k is the ra- 
tio of the string linear tension to the friction coefficient; 
(£(0)£(i)) = 2n8 (t) , and k is the thermal (molecular) 
diffusion coefficient. The solution of ([!]) is given by 



Pit) 



-kl 



W(t) 



P(0) 



e kt'W- 1 (t')W) 



I) 



(2) 



with W(t) is the time-ordered exponential 
- 1 i 

Texp J dt'a(t') . The temporal evolution of W is de- 
Lo J 

scribed by two Lyapunov exponents at times greater 
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than the inverse of the exponents, while all other ("an- 
gular") degrees of freedom are frozen |l8| . Respec- 
tively, ir[WW+] grows as exp[2 Xt] and one gets the 
following estimation for the length of the polymer, 
p 2 {t) ~ max{l,exp[2(A- k)t]}K/\X- k\. Here the 
long-time statistics of the leading Lyapunov exponent 
(according to a matrix version of the large deviation 
theory is described by the weight 



dAexp [-tS(A)] 



(3) 



The entropy function 5(A) is concave with a zero mini- 
mum at A = A. At the largest times, the moments of p are 
governed by a saddle point value, A = A„ , 5' (A„) = n, 



(P(t) n ) 
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This result shows, particularly, that only the lowest mo- 
ments (and only if k > A ) are steady. However small the 
A (in comparison with k ) may be it can not prevent the 
higher moments from an asymptotic in time growth, the 
steady PDF of p decays algebraically [^(J with a large-p- 
cutoff growing in time. The threshold for the n-th mo- 
ment convergence is given by 



5'(A„)=n, k n = X n - 5(A„) /n. 



(5) 



k n as a function of n has a minimum, A , at n = . The 
n = 2 value of the threshold is associated with the con- 
vergence of the dilute polymer solution contribution into 
the overall strain rate, and therefore can be identified as 
the onset of the drag reduction (20| . The dependence of 
the onset value on the number of elements in a long poly- 
mer chain (i.e. on square equilibrium radius of gyration 
of the polymer) will be discussed later in the letter. A 
problem, however, is that the linear theory does not al- 
low to describe statistics of polymer (and, particularly, 
to find the steady length of the stretched polymer) above 
the drag reduction onset. 

Account for anharmonic elasticity, we are switching 
our attention to, solves the latter problem. The anhar- 
monic generalization of (jl|) to consider is 



dt P 



(t)pf = -9;c/(p)+f(i), 



(6) 



where U(p) is an arbitrary bounded potential. Isotropy 
dictates that the distribution function of p a should de- 
pend on the polymer length only. For the dynami- 
cally relevant times, t ^> 1/A, the problem reduces to a 
study of (|^) projected on the major stretching direction. 
First, consider the case of the Gaussian entropy function, 
5(A) — (A — A) 2 /2A, for which the reduction from the 
stochastic differential equation (||) to the Fokker-Planck 
equation for the distribution function of the separation 
p, V, is especially simple 



dt + ^d p pd p p - Xd p p + 2nd 2 p + d p U\p) }V = 0. (7) 



The stationary solution of ((7J) is 



V(p) = Cp^exp 



- / dr 



, U'(r') - Ar' 
2k + Ar' 2 /2 



(8) 



where C is the normalization coefficient. One observes 
that when the molecular diffusivity is small it enters only 
as an ultraviolet cutoff, = ^/k/A in the large scale 
(p » rj ) behavior. In the forthcoming calculations the 
diffusion will be considered small. We will not be count- 
ing for effect of diffusion explicitly and yet introducing 
the Td cutoff if appropriate. The general case of a non- 
Gaussian entropy function does not lead to a simple dif- 
ferential equation for the PDF. However, it may always 
be described by the path integral 



p(0)=r 

V{r) = J VX{t)Vp(t)Vp(t)exp 
C = - P [p-Xp + U' (p)] + 5(A) 



Ldt 



, (9) 
(10) 



The effective action in (|9|) is large as the time of evolu- 
tion in the strong flow is. Therefore, the semi-classical 
saddle-point approximation should work. The saddle sys- 
tem of equations, supplied by the final condition (at the 
observation time) p(0) = r , is 



Psp ~ ^spPsp + U (p S p) — 0, 
~Psp — X S pP S p + p S pU" (p S p) = 0, 
PspPsp + 5' (A sp ) = 0. 



(11) 



Notice also that this system of equations is consistent 
with the conservation of energy 



E 



-p sp [-X S pp S p + U' (psp)] + 5 (A sp ) = const. (12) 



Because one is looking for finite action during an infinite 
time (beginning at t = — oo ) evolution, the choice E = 
is mandatory. The saddle point PDF is given by 



V sp (r) ~ exp 



dpsp 

Psp 



S' (X sp ) 



(13) 



with A sp and 



related to each other locally through 
U'(p S p)/ p S p — X sp - 5 (A sp _) /S' (Asp) . For the Gaussian 
5(A), one derives p sp - G + Xp sp -a - U' ipsp-a) = and 



Vsp-air) ~ exp 



2 
A 



% (F'{r>) - Ar') 



(14) 



Comparing the semi-classical expression ( |l4| ) with the ex- 
act one (pi), one concludes that the algebraic 1/r factor is 
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gained through the fluctuations calculated on the top of 
the saddle solution. (Notice also, that the saddle-point 
solution becomes exact after changing to log lnfr/rj vari- 
ables). To find the steady length of the polymer, one 
should also perform an additional variation of the effec- 
tive action in ( |l4| ) with respect to r. The result is given 
by the time-independent solution of (O) 



X Pat = U' (p st ) . (15) 
The last equation is valid for an arbitrary form of S(X). 

Consider now the case of a long polymer modeled by 
N anharmonic springs connected in a chain 



dt 



a a/3 (3 
Pi =° Pi- 



[d a U(p i;i ^)+d a U(p i;i+1 )]H?, (16) 



with pi-k = pi — pk , and the boundary condition, 
Pifl = Pn+i-.n = , imposed at the fictitious end points 
of the polymer, added for convenience. As the strong flow 
regime when the chain is extremely stretched (so that 
anharmonicity becomes more important while off-linear 
conformations of the polymer are less relevant) one can 
1) take into account only linear (straight-line) conforma- 
tions of the polymer; 2) neglect diffusion; 3) replace the 
projected (on the major stretching direction) a by the 
leading Lyapunov exponent. The continuous version of 
the scalar projection of (|l^) (to be considered for N 3> 1 , 
the major case of interest) is 



X(t) Pn +d n U' (d nP ), 



(17) 



supplemented by d n p (0) =d n p (N) = . The saddle- 
point method explained above in detail for the single- 
spring model is also valid (it is even enhanced as iV ap- 
pears to be an additional saddle parameter) for analyzing 
the following system of equations 



" X sp p sp + d n U' (dnPsp) = 0, 
-p sp - X S pP S p + d 2 p sp U" {d n p sp ) = 0, 

0. 



A' 



(18) 



/ PspPspdn + S' (A sp ) 
o 



zero energy condition 

-X a pp S p + d n U' {d n psp)} 



is 



The 

N 

J dnp S p 
o 

point system of equations supplied by the final condition 
p„(t = 0) = r(n) has a unique solution. The functional 



consistent with (|18| ) 

S(X sp ) . The saddle 



J — c 



i Psp Psp 



dt 



If. 



PDF is given by, V {r (n)} ~ exp 

however, one is interested in the PDF of the final point 
(polymer length) only, variation with respect to all the 
r(n) except the end of the polymer, r(N) , should be 
taken. A very symmetric solution of the enhanced sad- 
dle point system is available. The hint is that not only 
the global energy vanishes but also the energy density 
(for each and every value of n) does. This leads to 



, d n U'(d n p S p) _ S(X sp ) _ 

Psp /, n i 

s p 7, ~~ ~q77\ — V ~~ ' \ ' 

Psp >J (Asp) Psp 

The first equation in ( |l9| ) can be integrated 



I 



U" (x)dx 



n + f.V - 2n)V[—-n) = 



Jlfl* yU"(y)dy 



where r = d n p sp , and 0(x) is the step function. There- 
fore, the maximal tension on the string r*(t) and the 
length of the polymer, p*(t) — p sp (t;N), are related to 
each other through 



N* X sp - 



J .IT 



U"(x)dx 



pi 



S(Xsp) _ ^ 

S ' {Xsp) sfJl* yU"{y)dy 

S(X sp ) 



P S'(X sp ) 



xU"{x)dx. 



(20) 
(21) 



Finally, according to the second part of (19), applied 
to n = N , the dynamics of all the time-dependent 
fields ( X sp (t), p*{t) , T*(i) ) is unambiguously fixed by 
d\n\p*{t)]/dt = S(X sp )/S'(X sp ) , p*(0) = R . The prob- 
ability to observing a polymer of length R is given by ( |l3| ) 
with p sp replaced by , and A sp related to each other 
through ( pp| , ^l|) . The PDF has a pronounced maximum 
at some value, i?» , which is the steady state polymer 
length. i?» is described by (20 2]]) taken at A sp — * A 
(then p* R* ). 

The algebraic system of equations ( p0pl| ), with the re- 
sult substituted in the — > p sp version of (|T^), solves 
unambiguously the problem of finding the probability for 
the TV-chain polymer characterized by the inter-bead in- 
teraction U(x) to be stretched to the length R in the 
flow described by <S*(A). For the Gaussian entropy func- 
tion S (A) and the algebraic potential U(x) 
derives 



qx 



V(R) 



R* — 



2A/A 



exp 



2[X/ q ] 1/a 


( r 


4+3q 
4+2a 




(l + a) 1/a 


V r 


l+Q 
2+Q 


V2tt 




2+a 



(22) 
(23) 



Notice also, that (23) holds for an arbitrary form 
of 5(A). For the linear harmonic chain, U = 
kx 2 , and the Gaussian entropy function, one gets, 
A In [P{R)\ I (2 \n[R/r d ]) = X — ir 2 k/N 2 , i.e. for N 
greater then ir \fxjk the pure harmonic consideration has 
no sense as the resulted pdf is not normalizable. In other 
words, a transition to the nonlinear regime is taking place 
at a very small (ir 2 k/N 2 ) value of A. 
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Knowledge of P(R) allows to calculate any correlation 
functions of R. I will briefly discuss two examples of in- 
terest. First of all, as it was already mentioned above, 
the kinetic theory estimate for the polymer contribution 
into strain tensor is proportional to (R 2 ) = J R 2 P(R)dR 
ppf . Therefore, (|2^,^l|,0) applied to the case of a com- 
plex potential U{x) = kxr+qx 2+a , will describe the "first 
order phase transition" (jump of the (R 2 ) value) from the 
linear theory estimate, ~ r 2 ,, valid at A + 2A < Tr 2 k/N 2 , 
upto ~ R 2 , correspondent to the nonlinear case. Sec- 
ond example explains how to recalculate the probability, 
Pdcg, of the polymer to degrade (to break in two pieces) 
under the action of turbulence. Indeed, assume that the 
polymer string breaks if the maximal (along the polymer) 
tension exceeds the critical value, r c . Then, Pd eg is given 
by the integration of P{R) from R c upto oo, where R c 
solves the system ( po||2l] ) with p* and r* replaced by R c 
and r c respectively. 

Let me briefly mention other passive problems of inter- 
est which remained beyond the scope of the present let- 
ter. First of all, it is known that in the case of fully devel- 
oped turbulence, the degenerate (or very close to degen- 
erate) sheet-like flow configurations (associated with the 
so-called Vieillefosse tail of the velocity gradient PDF) 
can be favored, especially in the viscous range [El) . Poly- 
mer stretching in the degenerate case requires a separate 
investigation accounting for off-straight-line but planar 
conformations of the polymer. Notice also, that tur- 
bulent stretching of a flexible but inextensible polymer 
j| (the one which is considered to be a good model of 
B-DNA p2[|) is yet another problem which requires an 
accurate account for the off-straight-line conformations. 
Second, excluded volume effect along with hydrodynamic 
self-interaction of the polymer, known to be of a great im- 
portance in the equilibrium theory p3| , p4f , may also play 
a significant role in the random flow case. Third, com- 
pressibility of the flow (responsible for the sign change 
of A ^5|) may cause new effects interesting enough to 
study (the saddle point approach will not work if A is 
negative). Finally, steady statistics of comb-, star-like 
and randomly branched polymers, dendrimers or poly- 
meric membranes, placed in a chaotic flow can also be 
studied within the passive approach demonstrated above 
for the linear polymer. 

If concentration of polymer solution is small enough 
the passive approach developed above is valid. However, 
the elastic contribution into the strain tensor grows with 
the concentration. Once the elastic contribution into the 
stress tensor becomes of the order of the kinetic one the 
problem cannot be treated as passive. The steady length 
of the polymer may be still much smaller than the vis- 
cous scale. The scale separation, in principle, allows to 
generalize the passive approach of the present work and 
to describe various chaotic effects of the non-Newtonian 



hydrodynamics (see pq ] for an attempt to proceed in the 
active direction). 
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